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ABSTRACT 

The ability of a three-dimensional Navier-Stokes method, 
PAB3D, to predict nozzle afterbody flow at high Reynolds 
number was assessed. Predicted surface pressure coefficient 
distributions and integrated afterbody drag are compared with 
experimental data obtained from the NASA-Langley 0.3 m 
Transonic Cryogenic Tunnel. Predicted afterbody surface 
pressures matched experimental data fairly closely. The 
change in the pressure coefficient distribution with Reynolds 
number was slightly over-predicted. Integrated afterbody drag 
was typically high compared to the experimental data. The 
change in afterbody pressure drag with Reynolds number was 
fairly small. The predicted point of flow separation on the 
nozzle was slightly downstream of that observed from oil- 
flow data at low Reynolds numbers and had a very slight 
Reynolds number dependence, moving slightly further down- 
stream as Reynolds number increased. 

INTRODUCTION 

Potential differences between the data obtained through 
typically low Reynolds number scale model wind-tunnel test- 
ing and that of full scale has historically been a concern and 
the present advanced subsonic and high-speed research pro- 
grams have renewed this attention. In particular, concern 
about potential Reynolds number effects on drag has been 
focused in the area of propulsion-airframe integration (PAI) 
testing, where high afterbody slopes and long boundary-layer 
runs result in large viscous effects on the nozzle boattail. 
Additionally, extensive use of computational fluid dynamics 
(CFD), from incompressible panel methods through advanced 
Navier-Stokes methods, for preliminary design studies in 
these same transport programs, with no direct evidence of 
their ability to predict correctly the effects due to Reynolds 
number, have resulted in similar concerns. 
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Previous to the development of cryogenic test techniques 
for achieving high Reynolds numbers in wind tunnel test fa- 
cilities little fundamental research data had been available 
for the evaluation of any CFD methods for their ability to 
predict Reynolds number effects in PAI. However, a number 
of years ago, during the developmental phase of cryogenic 
testing techniques at the NASA Langley Research Center; 
two sets of simple axisymmetric nacelle models were built 
and tested in what was then known as the 1/3 m Pilot Tran- 
sonic Cryogenic Tunnel (now the 0.3 m Transonic Cryogenic 
Tunnel); thereby obtaining some of the first test data for a 
number of nozzle-boattail configurations over a wide 
Reynolds number range, refs. 1-4. The 8 inch long (forebody) 
geometry of most of these models was, by design, a scaled 
down version of a 48 inch long model that had been tested in 
the Langley 16-Foot Transonic Tunnel, thus providing addi- 
tional comparisons of quality and repeatability of wind tun- 
nel data between facilities, ref. 1 . 

The multiblock three-dimensional Navier-Stokes method 
PAB3D, refs. 5, 6, and 7, was used for all the calculations 
and utilized the algebraic Reynolds stress model of Shih, Zhu, 
and Lumley for turbulence simulation, refs. 8 and 9. Grid 
sensitivity was evaluated for afterbody pressure coefficients 
and integrated pressure drag using the performance package, 
ref. 10. Afterbody pressure coefficients and drag were com- 
pared with experimental data. 

NOMENCLATURE 

Amax maximum body cross-sectional area, 0.78540 in^ 
Cp> pressure drag coefficient, F/qooA max 

Cp pressure coefficient, (p - pooVqoo 

d m body maximum diameter, 1 .0 in. 

L length of body from nose to nozzle connect station, 

8 or 16 in. 

1 length of nozzle, in. 

Moo free-stream Mach number 

p static pressure 

Poo free-stream static pressure 

qoo free-stream dynamic pressure 



N R e Reynolds number based on body length 

Uqo free-stream velocity 

u axial velocity 

x axial distance 

y lateral distance from model centerline 

z vertical or radial distance from model centerline 

a flow angle-of-attack, degrees 

4> angular location of pressure orifices, degrees 

Subscripts 

P nozzle boattail component contribution 

CL centerline 

exp experiment 

t total conditions 

DISCUSSION 

Experimental Model and Data 

The model used for this study was one of six models 
that were built for the original Reynolds number study, refs. 1 
and 2. Four models with differing boattail geometry with a 
length of 8 inches from the nose to the start of the boattail 
(characteristic length) and two models with a length of 16 
inches were constructed. The boattail geometries had circu- 
lar arc (2), circular arc -conic, or contoured cross-sections. 
This investigation utilized the circular arc with a length-to- 
maximum-diameter ratio (fineness ratio - l/d m ) of 0.8 
boattail. Figure 1 is a photograph of the model mounted in 
the pilot tunnel. The nose of the model was a 28° cone, 1 .7956 
inches long fairing to the cylindrical body via a 1.3615 inch 
radius circular arc centered 2.125 inches downstream of the 
model leading edge and 0.8615 inches below the model 
centerline. The circular arc fairing is tangent at its endpoints 
to the conical nose (1.7956 inches from the nose) and cylin- 
drical body (2.125 inches from the nose). The model was 
sting mounted with the diameter of the sting being equal to 
the model base diameter. The length of the constant diam- 
eter portion of the sting (6.70 inches measured from the nozzle 
connect station) was such that, based on the work of Cahn, 
ref. 1 1 , there should be no effect of the sting flare down- 
stream of the nozzle trailing edge on the boattail pressure 
distributions. The following is a tabulation of the non-dimen- 
sional orifice locations. 


x/d m for L/Dm = 8 at 

-©- 

ti 

o 

o 

0 = 120° 

$ = 240° 

-0.2771 

-0.2761 

-0.2850 

-0.0256 

-0.0731 

-0.0700 

0.0770 

0.0256 

0.0345 

0.1765 

0.1287 

0.1270 

0.2750 

0.2257 

0.2260 

0.3679 

0.3240 

0.3279 

0.4675 

0.4180 

0.4200 

0.5749 

0.5166 

0.5220 

0.6698 

0.6165 

0.6376 

0.7746 

0.7280 

0.7400 

x/d m for L/Dm = 1 6 at 

-0.4491 

-0.4660 

-0.4561 

-0.1637 

-0.2201 

-0.1552 

-0.0600 

-0.1281 

-0.0590 

0.0337 

-0.0260 

0.0390 

0.1268 

0.0744 

0.1342 

0.2279 

0.1729 

0.2713 

0.3210 

0.2696 

0.3718 

0.4199 

0.3679 

0.4680 

0.5231 

0.4640 

0.5749 

0.6279 

0.6758 

0.7304 


The models were constructed of cast aluminum with 
stainless-steel pressure tubes cast as an integral part of the 
model. The models were instrumented with 30 pressure ori- 
fices in three rows of 10 orifices each. The 1 inch diameter 
of the model physically precluded the placement of all 30 
orifices along the same row. The tunnel has an octagonal test 
section with slots at the corners of the octagon and is essen- 
tially a scale model of the Langley 16-Foot Transonic Tun- 
nel test section, ref. 12. The test medium for the cryogenic 
tunnel was nitrogen. 

The total temperature and pressure of the pilot tunnel 
could be independently controlled. The experiment was con- 
ducted over a range of temperatures from approximately 1 1 7 
K to 308 K and pressures from 1 to 5 times the standard 
atmospheric level. Any number of free stream total tempera- 
tures or pressures can result in identical settings of Reynolds 
number. From experiments conducted at several conditions 
which resulted in the same Reynolds number surface pres- 
sure coefficients and nozzle boattail drag were shown to be 
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similar regardless of the temper ature/pressure combinations 
that created equivalent Reynolds numbers, ref. 4. Though 
data were obtained over the range of Mach number from 0,6 
to 0.9, only the M = 0.6 data are compared with the CFD of 
this paper. The following is a table of conditions for experi- 
mental data obtained at M = 0.6 for the L/d m = 8 model. One 
atmosphere is defined as 0.101325 MPa. 


T t ,K 

p t , atm 

N Re xicr 6 

117 

5.00 

43 

117 

4.00 

34 

117 

3.00 

26 

117 

2.50 

21.5 

117 

2.00 

17 

117 

1.50 

13 

117 

1.30 

1 1.5 

308 

5.00 

11.5 

308 

3.80 

8.5 

308 

3.14 

7.0 

308 

2.50 

5.5 

308 

1.25 

3.0 


Computational Procedure 

Governing Equations . The code used was the general 3-D 
Navier-Stokes method PAB3D-Version 13G. This code has 
several computational schemes and different turbulence mod- 
els that can be utilized, as described in more detail in refer- 
ences 6 and 8. The governing equations are the Reynolds- 
averaged simplified Navier-Stokes equations (RANS) 
obtained by neglecting all streamwise derivatives of the vis- 
cous terms. The resulting equations are written in general- 
ized coordinates and conservation form. The implementa- 
tion of the full three dimensional viscous stresses are reduced 
to thin-layer viscous assumptions, although full Navier-Stokes 
simulation is an option. The diffusion terms are centrally 
differenced and the inviscid flux terms are upwind 
differenced. Two finite volume flux-splitting schemes are 
used to construct the convective flux terms. The Roe upwind 
scheme with third order accuracy is used in evaluating the 
explicit part of the governing equations and the van Leer 
scheme is used to construct the implicit operator. 

The user will typically utilize the Roe scheme procedure 
to sweep streamwise through the computational domain and 
the van Leer scheme for the solution of the cross-plane (i.e M 
i - constant) of a three dimensional problem. A single-cell 


wide, two-dimensional mesh defined with the i direction of 
the grid oriented in the conventional streamwise direction will 
typically converge slower using the Roe relaxation solution 
scheme compared to solving the equivalent problem with the 
van Leer scheme. Therefore the i and j directions of the 2-D 
mesh are swapped allowing the entire flowfield to be solved 
implicitly with each iteration. The explicit sweep is not used 
since only one cell exists in the i direction. The implicit scheme 
has the potential of a much higher rate of convergence and can 
provide a solution using less computational time. 

Several near-wall models and compressibility corrections 
are available to be used with several formulations of linear 
two-equation k-£ turbulence model, (e.g. Standard, Jones & 
Launder, Yang & Shih). The k-£ turbulence model equations 
are uncoupled from the RANS equations and can be solved 
with a different time step than that of the principle flow solu- 
tion. Version 13G of the PAB3D code used in this study has 
options for several algebraic Reynolds stress (ASM) turbu- 
lence simulations. The standard model coefficients of the k- 
£ equations were used as the basis for all the linear and non- 
linear turbulent simulations, ref. 13. Additionally, it is known 
that the eddy viscosity models produce inaccurate normal 
Reynolds stresses. Flat plate flow, as well as other more com- 
plex aerodynamic flows, are anisotropic. Successful imple- 
mentation of the algebraic Reynolds stress models required 
the solution methodology for turbulent production term of 
the underlying linear turbulence calculations to be modified. 
Turbulence production depends on high order derivatives of 
the turbulent Reynolds stresses. Proper representation of the 
stresses should be provided by face centered values, rather 
than the cell centered values. Previous attempts to imple- 
ment non-linear turbulence models in the context of a cell 
centered eddy viscosity model worked only for 2-D prob- 
lems and was unable to resolve 3-D flows. The two non- 
linear turbulence models used in this paper are the theories 
by Shih, Zhu & Lumley, ref. 9, and by Gatski & Speziale, 
ref. 14. All flow solutions were developed with PAB3D- 
V 1 3G using the standard formulation of linear k-£ as the ba- 
sis for the algebraic Reynolds stress calculation. The damp- 
ing function of Launder & Sharma, ref. 1 5, was used to control 
the near- wall behavior of k and £ was set to 2o((9Vk)/5n)^ at 
solid surfaces. Turbulent flow solutions using either the ASM 
or two-equation linear k-£ model requires 23 words per grid 
point. The code speed was 43 micro-seconds per grid point 
(Cray 2 time) solving ASM turbulent flow simulations. 

The conservative patch interface package of Pao & 
Abdol-Hamid, ref. 7, enables the code to properly transmit 
information between mis-matched block interfaces. Integer- 




to-one interfaces are considered a subset of the arbitrary block 
interface and do not need to be specified as such to the patch- 
ing code. The patching program writes a connectivity data- 
base as a preprocessor prior to execution of the flow solver. 
Each entry to the patch database contains cell face areas and 
indices relating that cell with all other cells that will share 
momentum flux information. The database information is 
automatically re-allocated internal to the code during mesh 
sequencing. As a result, each block can be sequenced at dif- 
ferent levels and the correct interface information is main- 
tained at the cell level. However, it is important to note that 
features in the flow developed on one side of an interface 
should not be obliterated on the other side due to too severe 
a grid density mis-match. 

Third-order continuity in transmitting the fluxes across 
block boundaries is maintained by the code; lower order con- 
tinuity may be specified by the user if required. Equal cell 
size spacing on either side of an interface in directions nor- 
mal to the interface must still be maintained regardless of the 
mesh sequencing level of the block. 

Boundary Conditions . For this study, solid walls were treated 
as no-slip adiabatic surfaces. The solid wall boundary condi- 
tion was satisfied by setting the momentum flux of the solid 
wall cell face to zero. The boundary conditions used for the 
internal nozzle flow path were fixed total pressure, total tem- 
perature and flow angle at the plenum block inflow face. A 
boundary condition for the Riemann invariants along the 
characteristics was specified for the external freestream in- 
flow face and the lateral freestream outer boundary of the 
flow domain. An extrapolation boundary condition was ap- 
plied on the downstream outflow face where both the free 
stream and the nozzle plume exit the computational domain. 
The axisymmetric flow assumption for the single-cell grids 
was implemented by placing flow symmetry conditions to 
the lateral side boundaries of the computational domain. 

Solution Process. Total afterbody drag, nozzle pressure drag, 
and solution residual were used to determine the solution 
progress at the coarse (144), medium (122), and fine (111) 
grid levels of the axisymmetric afterbody. The 144 abbre- 
viation means the number of i-cells, j-cells, and k-cells were, 
internal to the code, reduced by 1,4, and 4 respectively to 
develop a coarse grid solution. Afterbody drag variance of 
less than 0.50 percent for several hundred iterations was 
achieved for all test cases. 

Computational Grid Setup . The axisymmetric afterbody grid 
utilized H-O type mesh topology with all block dimensions 


that were sequencable by 4. The mesh was gridded with a 
single cell wide 5 degree wedge grid with the streamwise 
flow direction oriented along the j index to utilize the im- 
plicit flow solver in the code for faster solution convergence. 
Figures 2 through 4 are overall and detail views of the cell 
centered grid. The body was described using 100 cells ex- 
tending from the leading edge of the nose to the nozzle con- 
nect station. There were 80 cells extending from the nozzle 
connect station to the nozzle boattail trailing edge. 

The first cell height typically decreased with increasing 
Reynolds number of the freestream. The following schedule 
was used for the freestream conditions of M = 0.6 and 
T t = 560° R and the model reference length of L/d m = 8. 
These data are also plotted in figure 5 and are compared with 
a 1-D isentropic flat plate curve of the variation of the first 
cell height with Reynolds number for a flat plate (solid line) 
for a y + = 1.0. Therefore, the flow in the first cell center 
above the surface using a grid generated with guidance from 
flat plate theory would then be y + = 0.5. 


N Re x 10- 6 

PlO* P si 

h j , in x 10 6 

y + 

11 

73. 

20 

0.50 

25.6 

171. 

8 

0.27 

43.1 

287. 

2 

0.18 

100.0 

665. 

2 

0.39 


The first physical cell height for the different Reynolds 
number grids for this model are indicated by the solid circles. 
The open circle symbols are the y+ of the flow and show that 
the y + for converged solutions were generally between 0.5 
and 0.18. The low Reynolds number grid calculation was 
done at y + = 0.5. The ASM turbulence simulations have bet- 
ter convergence trends when at least 2 cells are below y + = 1 . 
Therefore, for the higher Reynolds number grids the first cell 
spacing was slightly below the flat plate predictions. The 
range of acceptability is generally between a y + of 0.2 and 
0.5, hence the two higher Reynolds number grids actually 
had the same physical spacing but the y + was still inside the 
acceptable range. 

Results 

Figures 6 through 10 are the sensitivity of the CFD 
method to grid density with Reynolds number. Figure 1 1 is 
the change in predicted nozzle pressure coefficient distribu- 
tion with Reynolds number. Figure 1 2 shows the experimen- 
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tal trend in nozzle pressure coefficients with Reynolds num- 
ber. Figure 13 is a comparison of CFD with experimental 
data obtained in two different wind tunnels around the 
Reynolds number of 5 million. Figure 14 is a comparison of 
predicted and experimental nozzle pressure drag with 
Reynolds number and figure 15 is a comparison of two ana- 
lytic techniques for predicting skin friction drag on the nozzle 
boattail with Reynolds number. Figure 16 is the predicted 
trend of the point flow separation on the nozzle with Reynolds 
number. Figure 17 is a photograph of an oil-flow at M = 0.85 
and figure 18 is an analysis of integrated pressure drag error 
build-up between experimental data and CFD. 

In figures 6 through 10 repeat experimental data points 
were plotted along with solutions at three grid levels (coarse 
( 1 44), medium (122), and fine ( 1 1 1 )) to evaluate the experi- 
mental data scatter and CFD grid sensitivity, respectively. 
The coarse grid computation was consistently low in the re- 
covery pressure in the region of separated flow downstream 
of x/d m = 0.55 for the complete Reynolds number range. 
Grid convergence was closely attained as the medium and 
fine grids were very similar, but not identical, at each of the 
Reynolds numbers. The level of flow expansion around the 
boattail shoulder was predicted. The pressure recovery was 
generally underpredicted at Reynolds numbers below 40 
million, while the pressure level was generally matched above 
40 million. 

Figures 1 1 and 12 show the predicted change in nozzle 
surface pressure coefficients with Reynolds number 
computationally and experimentally, respectively. The ex- 
perimental data curves are smoothed, averaged data over the 
three angular rows of orifices. In general, the CFD shows 
lower peak pressure coefficients at the flow expansion peak 
near x/d m =0.15 with increasing Reynolds number. Also 
higher recovery pressures in the region of separated flow 
downstream of approximately x/d m = 0.6 are predicted with 
increasing Reynolds number. These trends were qualitatively 
the same with the experimental data, though the CFD over- 
predicted the change in recovery pressure in the separated 
flow region with Reynolds number. 

A comparison of CFD (computed for the 5.7 million 
Reynolds number of the cryogenic tunnel data) with experi- 
mental data obtained from both the 16-foot Transonic Tun- 
nel and the cryogenic Tunnel and two different wind tunnel 
models, ref. 16, is shown in figure 13. The differences in the 
flow due to the Reynolds number variation 3.6 million to 
6.9 million should not effect the validity of the comparison. 
There are data for the 1 6 inch model from both tunnels, refs. 


1 and 16, and data for the 8 inch model from just the 0.3 m 
tunnel, ref. 1 . In general there does not appear to be a signifi- 
cant facility related difference in the pressure coefficient dis- 
tributions. As mentioned in the previous section concerning 
grid sensitivity at the different Reynolds numbers, the CFD 
solution captures the flow expansion around the boattail 
shoulder fairly well, but slightly under-predicts the pressure 
recovery in the separated flow region downstream of 
x/d m = 0.5 at this low a Reynolds number. 

Figure 14 is a comparison of predicted nozzle boattail 
pressure drag coefficients with several experimental data sets. 
Like before, data from both the 0.3 m and 16-foot tunnels for 
several different models with the same nozzle boattail ge- 
ometry are presented. Additionally, data from the 48 inch 
model with the same boattail geometry, ref. 17, is included. 
In general, the predicted pressure drag was higher than the 
data obtained from the 0.3 m facility, ref. 1 . A better com- 
parison is obtained with the 16-foot data which are the open 
triangle, diamond, (ref. 16) and x symbols, ref. 17, at 11 mil- 
lion Reynolds number and below. The open diamond in fig- 
ure 14 is the integrated pressure drag obtained from the open 
circle symbols of figure 13. The open square around 6 mil- 
lion Reynolds number in figure 14 is the integrated pressure 
drag obtained from the solid circle symbols of figure 13, 
which appear fairly interspersed with the open circles. Ref- 
erence 1 noted that very slight differences in pressure coeffi- 
cient distributions appear to induce significant shifts in the 
level of pressure drag. The issue of assigning areas for the 
pressure-area integration process will be dealt with in a sub- 
sequent section. Also included in Figure 14 is a comparison 
of predicted nozzle boattail pressure drag for the 105 million 
Reynolds number test condition utilizing two alternative tur- 
bulence models, a linear model and that of Gatski and 
Speziale. There was not a significant difference in the drag 
coefficients predicted by the three methods. 

A related issue when comparing CFD with experimen- 
tal data is the assessment of skin friction drag. Typically the 
experiment assigns a wetted area to a portion of the model 
and, with the assumptions of constant Mach number and 
measured free-stream conditions, a flat-plate equivalent skin 
friction coefficient is calculated. The computational method 
has available considerably more information concerning the 
local flow on the model and could theoretically make an ex- 
act calculation of the skin friction forces. Three significant 
factors are accounted for in the computational method. First 
the stream wise variation in the local viscous stress at the 
wall can be calculated. Second, the effect of flow separation 
occurring on the latter 40 percent of the boattail can be bet- 
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ter accounted. Lastly, due to the curvature of the boattail, the 
component forces of the skin friction can be separately ac- 
counted. Despite these differences, the calculated skin-fric- 
tion drag compared to the estimation from the wind tunnel 
investigation was roughly 5 boattail drag counts off at 5 mil- 
lion Reynolds number and 2 boattail drag counts off at 100 
million Reynolds number, shown in figure 15. 

Figure 16 is a comparison of predicted point of flow sepa- 
ration on the boattail compared to flow visualization data 
obtained in the 16-Foot tunnel. Flow visualization for this 
model were only available at approximately 12 million 
Reynolds number. The 8, 16, and 48 inch models were run 
computationally and are shown by the solid symbols. The 
point of flow separation experimentally was determined from 
picking points from photographs of oil-flow runs in the 16- 
Ft. tunnel. The point of flow separation computationally was 
the sign reversal of the local skin friction coefficient. Though 
no statement can be made concerning the experimental trend 
of the point of flow separation with Reynolds number for 
this model, computationally, a slight downstream shift in the 
separation point is observed with increasing Reynolds num- 
ber. At 12 million Reynolds number, the CFD is slightly 
downstream of the experimentally estimated point. Figure 

1 7 is a photograph of the oil flow over the solid plume model 
at M = 0.6. As can be observed by the mildly wandering dark 
line part way down the nozzle boattail, some three-dimen- 
sionality existed in this flow that was not in the computa- 
tional solution. The dark line is assumed to be the flow sepa- 
ration point and the variation in location was estimated and 
is denoted by the vertical spread in open triangle symbols in 
figure 16. 

As an attempt to clarify the difference between the inte- 
grated pressure drag on the boattail computationally com- 
pared with the experimental data when the pressure coeffi- 
cient distributions appear fairly similar, several methods of 
integrating the computational data were investigated. The 
CFD drag numbers plotted in figure 14 were obtained by 
using the areas and pressures obtained directly from the com- 
putational grid. Due to the availability of the experimental 
areas used for the pressure integration, tabulated below, in- 
terpolated CFD pressure data were obtained at the exact ex- 
perimental orifice locations and subsequently integrated us- 
ing the experimental areas. The open circle symbols in figure 

1 8 are the experimental data and the x symbols are the inter- 
polated CFD solutions at the same physical locations. The 
filled diamonds are the difference between the CFD and ex- 
perimental pressure coefficients multiplied by the axial area 
for that orifice location. Several locations have a large dif- 


ference between the CFD and experiment but a very small 
projected area so the combination makes little difference in 
the error. The flow further downstream, between x/d m = 0.2 
to 0.5, experiences a rapid recompression trend that, though 
there is only a very slight downstream shift in the predicted 
recompression, created a consistently positive increment be- 
tween the data and CFD. This difference multiplied by the 
increasing axial area contribution biased the predicted pres- 
sure drag. Overall, while there appeared to be a good com- 
parison between the pressure coefficient data, the integrated 
drag numbers can be substantially different. 


Areas for L/D m = 1 6 at 

<j) = 0° 

o 

O 

CM 

il 

4> = 240° 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.00051 

0.00051 

0.00051 

0.0 

0.00999 

0.01366 

0.02354 

0.01858 

0.01346 

0.03332 

0.01346 

0.01880 

0.04049 

0.01880 

0.02369 

0.05384 

0.02369 

0.05487 

0.05493 

0.05705 

0.10745 


A comparison of the integrated drags obtained from ex- 
periment, CFD, and CFD using the experimental areas for 
four values of Reynolds number can be found in the follow- 
ing table. As can be seen, the CFD consistently predicted 
highter drag than the experiment, no matter which area as- 
signment was used. 


N Re x 10- 6 

C D 

^exp 

c Dcfd 

c DcfD 

exp area 

25.6 

0.0123 

0.0214 

0.0218 

43.1 

0.0100 

0.0208 

0.0274 

67.0 

0.0118 

0.0184 

0.0194 

105.0 

0.0102 

0.0181 

0.0189 


However, it must be pointed out that these drag values 
are for an isolated nacelle and when applied to an actual air- 
craft the drag coefficients would be relatively lower. If the 
rule of thumb that when applied to a typical fighter the drag 
from an isolated nacelle should be divided by 20 is used the 
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next table results. With one exception, the drag predicted by 
CFD is within 5 airplane drag counts of that from the experi- 
ment, reasonably good agreement. 


N Re xlO- 3 * * 6 

^^exp 

c DcFD 

Cdcfd 

exp area 

25.6 

0.0006 

0.0011 

0.001 1 

43.1 

0.0005 

0.0010 

0.0014 

67.0 

0.0006 

0.0009 

0.0010 

105.0 

0.0005 

0.0009 

0.0009 


■SUMMARY 


4. Reubush, David E.: Experimental Investigation to 
Validate Use of Cryogenic Temperatures to Achieve High 
Reynolds Numbers in Boattail Pressure Testing. NASA TM 
X-3396, August, 1976. 

5. Abdol-Hamid, K. S.: A Multiblock/Multizone Code 
(PAB3D-v2) for the Three-Dimensional Navier-Stokes 
Equations: Preliminary Applications, 

NASA CR- 182032, October 1990. 

6. Abdol-Hamid, K. S.; Carlson, J. R.; and Lakshmanan, B.: 
Application of Navier-Stokes Code PAB3D to Attached and 
Separated Flows for Use With k-e Turbulence Model, NASA 
TP- 3480, 1994. 


The ability of a three-dimensional Navier-Stokes method, 
PAB3D, to predict nozzle afterbody flow at high Reynolds 
number was assessed. Predicted surface pressure coefficient 
distributions and integrated afterbody drag are compared with 
experimental data obtained primarily from the NASA-Lan- 
gley 0.3 m Transonic Cryogenic Tunnel with some data from 
the Langley 16-Foot Transonic Tunnel. Predicted afterbody 
surface pressures matched experimental data fairly closely. 
The change in the pressure coefficient distribution with 
Reynolds number was slightly over-predicted. Integrated 
afterbody drag was typically high compared to the experi- 
mental data. The change in afterbody pressure drag with 
Reynolds number was fairly small. The predicted point of 
flow separation on the nozzle was slightly downstream of 
that observed from oil-flow data at low Reynolds numbers 
and had a very slight Reynolds number dependence, moving 
slightly further downstream as Reynolds number increased. 
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Figure 3.- Cell centered grid for model. 



Figure L- Model installed in cryogenic tunnel. 




Figure 2.- Overall view of cell centered grid. Figure 5.- Variation of First physical cell height (yj) and y + 

for converged solutions with Reynolds number for freestream 
conditions of M = 0.6, T t = 560° R, and model with L/d m = 8 
compared with variation for a 1-D isentropic flat plate. 
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PAB3D, SZL, L = 1 6, coarse grid 

PAB3D, SZL, L = 1 6, medium grid 

— — ■ PAB3D, SZL, L = 16, fine grid 



Figure 6.- Comparison of pressure coefficient distributions 
from solutions at 3 levels of grid density for configuration 
with L/d m = 1 6 with experimental data for a Reynolds number 
of about 5.8 million. 


PAB3D, SZL, L = 8, coarse grid 

PAB3D, SZL, L = 8, medium grid 

— - PAB3D, SZL, L = 8, fine grid 



Figure 7.- Comparison of pressure coefficient distributions 
from solutions at 3 levels of grid density for configuration 
with L/d m = 8 with experimental data for a Reynolds number 
of about 1 1 million. 


PAB3D, SZL, L = 8, coarse grid 

PAB3D, SZL, L = 8, medium grid 

— PAB3D, SZL, L = 8, fine grid 



Figure 8.- Comparison of pressure coefficient distributions 
from solutions at 3 levels of grid density for configuration 
with L/d m = 8 with experimental data for a Reynolds number 
of about 43 million. 

PAB3D, SZL, L = 16, coarse grid 

PAB3D, SZL, L = 16, medium grid 

— — PAB3D, SZL, L s 16, fine grid 



Figure 9.- Comparison of pressure coefficient distributions 
from solutions at 3 levels of grid density for configuration 
with L/d m = 16 with experimental data for a Reynolds number 
of about 67 million. 
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PAB3D, SZL, L = 16, coarse grid 

PAB3D, SZL, L = 16, medium grid 

— »— PAB3D, SZL, L * 16, fine grid 
• L = 16, M =0.61, 105.7M, 4> = 0. 
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Figure 10.- Comparison of pressure coefficent distributions 
from solutions at three levels of grid density for configuration 
with L/d m = 16 with experimental data for a Reynolds number 
of about 105 million. 



Figure 11.- Comparison of computed pressure coefficient 
distributions at four Reynolds numbers. 



Figure 12.- Comparison of smoothed, averaged experimental 
pressure coefficient distributions at two Reynolds numbers. 


— « — PAB3D, SZL, L * 16, coarse grid 

• 0.3m, L = 1 6, M = 0.61 , 5.7M, 0 = 0. 

• 0.3m, L = 1 6, M = 0.61 , 5.7M, 0 = 120. 

• 0.3m, L = 16, M = 0.61 , 5.7M, $ = 240. 

O 16 ft, L = 16, M « 0.60, 3.6M, 0 = 0 
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Figure 13.- Comparison of computed pressure coefficient 
distribution with experimental distributions from both the 
cryogenic and the 16-Foot Transonic tunnels at Reynolds 
numbers from 3.6 million to 6.9 million. 
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0.8 r 


• PAB3D, -L = 8, SZL, fine grid 
■ PAB3D, -L =16, SZL, fine grid 
a PAB3D, -L = 16, G&S, fine grid 

♦ PAB3D, -L = 16, linear, fine grid 
O 0.3m cryo. data, L = 8 

□ 0.3m cryo. data, L = 16 
a 16 ft data, L = 8 
o 16 ft data, L= 16 
x 16 ft data, L = 48 



Figure 14.- Comparison of computed afterbody drag 
coefficient variation with Reynolds number with several sets 
of experimental data. 


c D,p,sf 


c D,psf,fp- L = 8 

• c D,(3sf,CFD’ L = 8 

^D.psf.fp' L ~ 16 
■ c D,(3sf,CFD' L = 16 
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Figure 15.- Comparison of computed afterbody skin friction 
drag coefficient variation with Reynolds number for both 
CFD and flat plate methods (used with experimental pressure 
drag coefficient data). 
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• PAB3D, SZL, L = 8, M = 0.6 
■ PAB3D, SZL, L = 16, M = 0.6 
a PAB3D, SZL, L = 48, M = 0.6 
-a- 16 ft, L = 48, M = 0.6, solid plume 
x 16 ft, L - 48, M = 0.6, NPR = 2 
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Figure 1 6.- Comparison of predicted point of flow separation 
with that obtained from flow visualization in the 16-Foot 
Transonic Tunnel. 



Figure 17.- Oil flow photograph from the 16-Foot Transonic 
Tunnel. 
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Figure 18.- Comparison of computed and experimental 
pressure coefficients at the same locations on the afterbody 
at a Reynolds number of 105 million and an indication of the 
contributions of the differences to drag coefficient differences. 
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